Ab-initio calculations of Many-Body effects in liquids: 
the electronic excitations of water 
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We present ab-initio calculations of the excited state properties of liquid water in the framework 
of Many-Body Green's function formalism. Snapshots taken from molecular dynamics simulations 
are used as input geometries to calculate electronic and optical spectra, and the results are averaged 
over the different configurations. The optical absorption spectra with the inclusion of excitonic 
effects are calculated by solving the Bethe-Salpeter equation. These calculations are made possible 
by exploiting the insensitivity of screening effects to a particular configuration. The resulting spectra 
are strongly modified by many-body effects, both concerning peak energies and lineshapes, and are 
in good agreement with experiments. 

PACS numbers: 78.20.C1, 78.40.Dw, 78.40.Pg, 71.35.-y 



The electronic structure of liquid water is still not 
fully elucidated, even though it is essential to under- 
stand the chemical and physical properties of many bio- 
chemical and industrial processes that occur in solution, 
where it is crucial to include the role of the solvent in 
the reactions. Water is also essential, both pure and 
as a solvent, for living organism survival and for bio- 
logical systems in general. For these reasons the study 
of the excited state properties of liquid water is funda- 
mental to advance in many research fields. However, 
in the last years, the theoretical studies of liquid water 
[1110,1111110,1111113 have mostly focused on its 
structure and ground state properties whereas less effort 
has been dedicated to its electronic structure and optical 
absorption spectrum. As a consequence, experimental 
data about excited states are poorly understood. One of 
the purposes of the present work is to solve these issues 
by carrying out ab-initio many-body calculations of the 
electronic structure and optical spectra of liquid water. 

Here we generalize the application of the Many Body 
Perturbation Theory (MBPT) [ll| to hquid systems, and 
present a calculation of the optical absorption spectrum 
of liquid water from first principles, including both self- 
energy effects and the electron-hole interaction. We show 
the occurrence of important excitonic effects, which are 
crucial for a good description and interpretation of ex- 
perimental data. 

The main problem concerning the study of liquid water 
relies on the fact that, in order to simulate a complex 
disordered system, one should use a very large unit cell, 
with a consequent prohibitive computational cost. In 
order to overcome this bottleneck, we used a smaller 
unit cell but exploiting several molecular dynamics 
(MD) snapshots of water as input geometries. Averaging 



the resulting electronic and optical spectra over many 
configurations should give a good approximation for the 
excited states properties of the real system. We used 20 
configurations of 17 water molecules in a cubic box with 
15 a.u. side and 8 special k-points. This is a quite low 
number of molecules [l(J; still, the computational effort 
of the many-body calculations for all the configurations 
would have been almost prohibitive. We will show in the 
following that one can however restrict the calculation 
of certain ingredients, in particular of the screening, 
to a single configuration. In this way the calculations 
become feasible with a reasonable effort. Already the 
results obtained for 17 molecules, averaged over several 
snapshots, well compare with experiments. 

The water configurations have been obtained by sam- 
pling every 2 ns a 40 ns long classical MD simulation 
trajectory. A TIP3P water model potential 12] has 
been used to represent the water molecules in our sim- 
ulation box. Equations of motion have been integrated 
numerically using a time-step of 1 fs. The MD run has 
been done in the NVT ensemble, where thermal equilib- 
rium at 298 K has been achieved applying a Nose-Hoover 
thermostat 13). Despite the small size of our system, the 
resulting radial distribution functions are in very good 
agreement with the experimental ones Il4| , as it is shown 
in Fig. n This confirms that we are using good input 
geometries for the excited state properties calculations. 

We then performed the electronic and spectroscopic 
calculations on three levels: we started with density 
functional theory (DFT) [Hi, [Hi to obtain the Kohn- 
Sham (KS) eigenvalues and eigenvectors. We hence con- 
structed the DFT independent-particle absorption spec- 
trum. As expected, it shows strong discrepancies with 
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FIG. 1: (Color online) Oxygen-oxygen (blue line) and 
oxygen-hydrogen (red line) radial distribution functions av- 
eraged over the 20 molecular dynamics snapshots compared 
to experimental data [l4| (circles and squares). 
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FIG. 2: Left panel; DFT HOMO and LUMO energies, aver- 
aged over the 8 k-points, for the 20 MD configurations. Right 
panel: schematic HOMO-LUMO gaps, calculated within DFT 
and GW, averaged over the 20 MD configurations. 



respect to experiment. These could not be cured by mov- 
ing to Time-Dependent (TD) DFT, in the adiabatic LDA 
(ALDA)[l3]. Therefore, we corrected the Kohn-Sham en- 
ergy levels using the Green's function perturbation ap- 
proach, with the exchange-correlation self-energy S cal- 
culated within the GW approximation E = iGW |l8l|.|llj 
(G is the one-particle Green's function and W = £~^v the 
screened Coulomb interaction). The Quasi-Particle (QP) 
energies (that is, the electronic 'band structure' of water), 
were calculated in first order perturbation theory: 

sT = e^J"" + A^r = e,r^+ < «|S - K^^l- > (1) 

where V^^^ is the exchange and correlation KS potential 
and the self-energy is evaluated at the quasiparticle en- 
ergy. As a last step, the optical absorption spectrum was 
calculated by solving the Bethe-Salpeter equation [ll|, 
so fully including local-field effects and the electron-hole 
interaction [ijj. 

DFT-KS results for the electronic and spectroscopic 
properties of water (obtained averaging over the 8 k 
points for each of the 20 configurations) are shown in 
Figs. [3 and ^. The configuration- averaged HOMO- 
LUMO gap turned out to be 4.85 eV, in good agreement 
with previous DFT calculations [l| but strongly under- 
estimating the experimental gap (8.7 ±0.5 eV [2i|), as 
expected in DFT calculations. Also optical absorption 
spectra at the DFT independent-particle level obtained 
for the 20 MD snapshots and, more important, their av- 
erage (Fig. |3|) do not compare in a satisfactory way with 
the experimental absorption spectrum J21| , |22| reported 
in the inset of Fig|Sl we can observe that the onset of the 
absorption is strongly underestimated in our calculation, 
the peak positions are red-shifted in comparison to the 
experiment, and the relative intensities of the first two 
absorption peaks are not well reproduced. 

We have therefore performed calculations of the opti- 
cal absorption spectrum within TDDFT |l7l |. In prin- 
ciple this represents an exact way to calculate optical 



spectra, but the quality of the results depends on the ap- 
proximation that is used to describe exchange-correlation 
effects. A widely used and computationally most effi- 
cient approximation is ALDA. Unfortunately our results, 
shown in Fig. 2Ji, shows no improvement with respect to 
the static DFT independent-particle result. Long-range 
and/or dynamical effects that are missing in the ALDA 
kernel should hence be important for the absorption spec- 
trum of water, and one has to resort to more elaborate 
(hence time-consuming) TDDFT approximations [23J , or 
to work in a different framework. We have adopted the 
latter choice in this work, since Many-Body Perturba- 
tion Theory offers a well established way to compute and 
analyze optical absorption spectra of extended systems. 

In order to correct the KS electronic gap and opti- 
cal spectra, one can calculate the GW corrections Ae^^ 
to the KS energies. This should be done for all the 20 
MD configuration, followed by an average. In particular 
the calculation of the screened Coulomb interaction for 
20 configurations constitutes however a true bottleneck. 
Instead, one can imagine that a change in configuration 
does not modify drastically the screening of an additional 
hole or electron. In fact, changing for example the posi- 
tion in space of an occupied orbital also moves the region 
of strong screening in the same direction. One can there- 
fore hope that GW corrections Ae^^ are quite stable 
with respect to the configuration. 

This actually turned out to be true in our calculations. 
In other words, although the configurations were differ- 
ent, and so were the KS eigenvalues, the difference be- 
tween DFT and GW, AGW = e^^ - e^-^^, was practi- 
cally constant going from one snapshot to another. This 
is shown in tabic Q for three configurations. We could, 
hence, use the same GW corrections for all the DFT con- 
figurations. 

The GW corrections increase the electronic HOMO- 
LUMO gap to 8.6 eV (see Fig. ^ "25|, well within the 
experimental range '20|. (Note that a larger gap, ~ 9.5 
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1.82 


3.86 



TABLE I: Test GW corrections to the HOMO and LUMO 
energy levels and to the electronic gap, for three different 
water configurations (E19, EOl, E13) Q^. Energies in eV. 




FIG. 3: Optical absorption spectra of the 20 configurations 
(thin lines) and their average (thick lines), a): DFT single 
particle calculation; b): with the inclusion of self-energy and 
electron-hole interaction effects. 



eV, is obtained in hexagonal ice |23. The disordered 
band-edge landscape leads indeed to a shrinkage of the 
niinimum gap with respect to the reference ordered sys- 
tem, ice in this case.). Since the calculated GW shift is 
almost constant for all the bands, GW optical spectra 
show lineshapes very similar to DFT ones, but shifted 
to higher energies, as shown in Fig. |31 The agreement 
with the experimental S2{uj) is hence not at all improved, 
since the relative intensities of the two absorption struc- 
tures (seen in experiments at 8.3 and 9.6 eV) are still not 
reproduced, and their position, from being red shifted 
in DFT calculations (Fig|31 solid line), are now strongly 
blue shifted in GW (FigEl dotted line). 

We have hence clear hints that single particle and sin- 
gle quasi-particle approaches are not sufficient to describe 
the optical properties of water, and that it is necessary 
to include the electron-hole interaction. To this end we 
had to solve the Bethe Salpeter equation, where electrons 
and holes interact through the screened Coulomb poten- 
tial W. These cumbersome calculations had to be done 
for all the 20 MD configurations, in order to obtain an 
average. 

Our results are shown in Fig. Ofor all the 20 MD con- 




FIG. 4: Optical absorption spectrum for one MD configu- 
ration, a) Time Dependent DFT within the Adiabatic LDA 
vs single particle Kohn-Sham DFT spectrum, b) BSE spec- 
trum calculated using the full screened electron-hole interac- 
tion (solid line) and using for W a constant dielectric function 
£"^ = 1/1.7 (dashed line). 



figurations and in Fig. [SJ where their average is reported 
(dashed line). Dramatic many-body effects are evident. 
Agreement with experiment both in energy peak posi- 
tions and onset, as well as in the relative intensities of 
the first two peaks, is significantly improved. The main 
remaining discrepancy is an overall redshift, that might 
be due to the fact that our GW calculations are not self- 
consistent but use DFT wavefunctions and energies. The 
first peak in the spectrum turns out to be a bound exci- 
ton with a binding energy of 2.4 eV and large oscillator 
strength. These are a consequence of the weak electronic 
dielectric screening of water (soo ~ 1.8). The second 
peak results from an excitonic enhancement of the os- 
cillator strength of interband transitions with respect to 
the single quasi-particle case. The binding energy of the 
lower exciton, although quite large, is smaller than the 
value Eb— 3.2 eV found in hexagonal ice 26] . One can in 
fact imagine that the mixing of electron-hole pairs of dif- 
ferent energy, which leads to the bound exciton, partially 
counterbalances the disorder effect on the quasiparticle 
gap [22] (namely the gap shrinking that is determined by 
a local increase of occupied state energies, and a local 
decrease of unoccupied state energies). Moreover, the 
higher density in the liquid phase with respect to ice 
(around 7%), may also play a role in reducing the ex- 
citon binding energy, by allowing a greater mobility of 
electrons and holes. 

It is of fundamental importance to explore also to 
which extent a detailed calculation of each single config- 
uration is mandatory. In fact, the effect of the electron- 
hole interaction does not depend significantly on de- 
tails of the screened Coulomb interaction W = e~^v 
[llj, but rather on macroscopic averages. Therefore, 
and as we have explicitly verified, the optical spectra 
of the various MD configurations can be calculated us- 
ing for the electron-hole interaction the same screening. 
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are well reproduced only if the electron-hole interaction 
is fully taken into account. ALDA calculations do not 
improve at all the DFT spectra. 
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FIG. 5: Calculated optical absorption spectra averaged over 
20 configurations. Solid line: single particle DFT spectrum; 
dotted line: spectrum obtained with the inclusion of the GW 
corrections; dashed line: spectrum obtained with the inclu- 
sion of the electron-hole interaction. Inset: experimental op- 
tical absorption spectrum of water (line: |2l|l : circles: |23|'). 



e~-^{q,uj)G,Gf, obtained for any of the independent MD 
snapshots. Moreover, test calculations done using a con- 
stant value for e^^ have shown that a good agreement can 
also be obtained using a value for e~^ = 1/1.7, (very near 
to the experimental macroscopic electronic value 1/1.8) 
as shown in Fig^. 

Our findings concerning the stability of GW and BSE 
calculations may have important consequences on future 
calculations for water (and maybe for other liquids), since 
the ab-initio evaluation of e~^{q,uj)G,G' is the real bot- 
tleneck of many-body calculations, and its evaluation 
for many different snapshots makes this approach too 
cumbersome for becoming a state of the art method. 
Hence, to be able to determine accurate optical spectra 
for many configurations, using the same dielectric con- 
stant to screen the electron-hole interaction, will enor- 
mously speed up the calculations. 

In conclusion, we have presented an ab-initio many- 
body calculation of the electronic structure and optical 
properties of liquid water obtained averaging the results 
of snapshots taken from molecular dynamics simulation. 
We have found that the GW corrections are almost inde- 
pendent of the particular snapshot considered and give an 
average HOMO-LUMO electronic gap of 8.6 eV, smaller 
than that of hexagonal ice [26j, and that the effect of 
the screened electron-hole interaction is the same for all 
the configurations. This result implies a huge reduction 
of the computational effort, and opens a pathway to low 
cost calculations on other disordered systems. 

Absorption spectra calculated with the inclusion of ex- 
citonic effects show important structures related to the 
electron-hole interaction which are essential to obtain a 
good description of experimental data. The onset po- 
sition and the relative intensities of the first two peaks 
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